Superconducting cosmic strings as Gamma Ray Burst engines 
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Cusps of superconducting strings can serve as GRB engines. A powerful beamed pulse of elec- 
tromagnetic radiation from a cusp produces a jet of accelerated particles, whose propagation is 
terminated by the shock responsible for GRB. A single free parameter, the string scale of symme- 
try breaking ~ 10^'* GeV, together with reasonable assumptions about the magnitude of cosmic 
magnetic fields and the fraction of volume that they occupy, explains the GRB rate, duration and 
fluence, as well as the observed ranges of these quantities. The wiggles on the string can drive the 
short-time structures of GRB. This model predicts that GRBs are accompanied by strong bursts of 
gravitational radiation which should be detectable by LIGO, VIRGO and LISA detectors. 



Models of gamma ray bursts (GRBs) face the prob- 
lem of explaining the tremendous energy released by the 
central engine In the case of isotropic emission, the 
total energy output should be as high as 4 x 10^** ergs. 
Strongly beamed emission is needed for all known engine 
models, such as mergers and hypernovae, but such ex- 
treme beaming is difficult to arrange. In this paper we 
show that emission of pulsed electromagnetic radiation 
from cusps of superconducting cosmic strings naturally 
solves this problem and explains the observational GRB 
data using only one engine parameter. 

Cosmic strings are linear defects that could be formed 
at a symmetry breaking phase transition in the early uni- 
verse 1^. Strings predicted in most grand unified mod- 
els respond to external electromagnetic fields as thin su- 
perconducting wires As they move through cosmic 
magnetic fields, such strings develop electric currents. 
Oscillating loops of superconducting string emit short 
bursts of highly beamed electromagnetic radiation and 
high-energy particles [^,^ . 

The idea that GRBs could be produced at cusps of 
superconducting strings was first suggested by Babul, 
Paczynski and Spergel (BPS) and further explored by 
Paczynski . They assumed that the bursts originate at 
very high redshifts (z - 100 - 1000), with GRB photons 
produced either directly or in electromagnetic cascades 
developing due to interaction with the microwave back- 
ground. This model requires the existence of a strong 
primordial magnetic field to generate the string currents. 

As it stands, the BPS model does not agree with 
observations. The observed GRB redshifts are in the 
range z ^ 3, and the observed duration of the bursts 
(lO^^s ^ T < 10*^5) is significantly longer than that pre- 
dicted by the model. On the theoretical side, our un- 
derstanding of cosmic string evolution and of the GRB 
generation in relativistic jets have considerably evolved 
since the BPS papers were written. Our goal in this pa- 
per is to revive the BPS idea taking stock of these recent 
advances. 



As in the BPS model we shall use the cusp of a super- 
conducting string as the central engine in GRB. It pro- 
vides the tremendous engine energy naturally beamed. 
Our main observation is that putting superconducting 
cusps in a different enviroment, the magnetized plasma 
at a relatively small redshift z, results in a different mech- 
anism of gamma radiation, which leads to a good agree- 
ment with GRB observational data. 

GRB radiation in our model arises as follows. Low- 
frequency electromagnetic radiation from a cusp loses its 
energy by accelerating particles of the plasma to very 
large Lorentz factors. Like the initial electromagnetic 
pulse, the particles are beamed and give rise to a hydro- 
dynamical flow in the surrounding gas, terminated by 
a shock, as in the standard fireball theory of GRB ||] 
(for a review see [|l[). We shall assume that cosmic mag- 
netic fields were generated at moderate redshifts {e.g., 
in young galaxies during the bright phase of their evolu- 
tion [^). The string symmetry breaking scale rj will be 
the only string parameter used in our calculations. With 
reasonable assumptions about the magnitude of cosmic 
magnetic fields and the fraction of volume in the Universe 
that they occupy, this parameter is sufficient to account 
for all main GRB observational quantities: the duration 
TGRB, the rate of events Ngrb, and the fluence S. 

We begin with the description of some properties of 
strings, which will be used below (for a review see j^]). 

A horizon-size volume at any time t contains a few long 
strings stretching across the volume and a large number 
of small closed loops. The typical length of a loop at cos- 
mological time t and the loop number density are given 
by 



I ~ at. 



ni{t) 



(1) 



The exact value of the parameter a in (|I|) is not known. 
We shall assume, following jl^], that a is determined 
by the gravitational backreaction, so that a kgGfi, 
where kg ~ 50 is a numerical coefRcient, G is Newton's 
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constant, ~ 77^ is the mass per unit length of string, 
and ry is the symmetry breaking scale of strings. 

The loops oscillate and lose their energy, mostly by 
gravitational radiation. For a loop of invariant length I 
[ pl| , the oscillation period is Ti ^ 1/2 and the lifetime is 
Ti ~ l/kgGfj,. 

An electric field E applied along a superconducting 
string generates an electric current. A superconducting 
loop of string oscillating in a magnetic field B acts as an 
ac generator and develops an ac current of amplitude 

Jo - e^Bl. (2) 

The local value of the current in the loop can be greatly 
enhanced in near-cusp regions where, for a short period 
of time, the string reaches a speed very close to the speed 
of light. Cusps tend to be formed a few times during each 
osciUation period. Near a cusp, the string gets contracted 
by a large factor, its rest energy being turned into kinetic 
energy. The density of charge carriers, and thus the cur- 
rent, are enhanced by the same factor. The contraction 
factor increases as one approaches the point of the cusp. 

The growth of electric current near the cusp due to 
string contraction is terminated at a critical value Jmax 
when the energy of charge carriers becomes comparable 
to that of the string itself, (J/e)^ ~ fi. This gives 

(erj/Jo). (3) 

Alternatively, the cusp development can be terminated 
by small-scale wiggles on the string ||l2| . If the wiggles 
contribute a fraction e <C 1 to the total energy of the 
string, then the maximum Lorentz factor is less than (|^), 
and is given by 'jmax ~ e~^^'^. The actual value of ^„iax 
is not important for most of the following discussion. 

Due to the large current, the cusp produces a pow- 
erful pulse of electromagnetic radiation. The total en- 
ergy of the pulse is given by |@,|| 5*^ 2kemJoJniaxl, 
where I ^ at is the length of the loop, and the coeffi- 
cient kern ^ 10 is taken from numerical calculations [Q. 
This radiation is emitted within a very narrow cone of 
opcnening angle 6min ~ 1 /jmax (for relativistic beaming 
in GRB sec fl^). The angular distribution of radiated 
energy at larger angles is given by Q 

d£em/dn ^ k,mJli/e\ (4) 

We shall adopt the following simple model of cosmic 
magnetic fields. We shall assume that magnetic fields 
were generated at some z ~ and then remained frozen 
in the extragalactic plasma, with 

= Bo(l + z)^ (5) 

where Bq is the characteristic field strength at the present 
time. For numerical estimates below we shall use zb ~ 4, 
Bq ^ 10^'' G, and assume that the fraction of volume of 



the universe occupied by magnetized plasma is Jb 0.1. 
We shall also assume that the universe is spatially flat 
and is dominated by non-relativistic matter. 

We shall now estimate the physical quantities charac- 
terizing GRBs powered by superconducting strings. For 
a GRB originating at redshift z and seen at angle 6 with 
respect to the string velocity at the cusp, we have from 
Eqs.(D-(D 

d£em/dn ^ k^me^ah^B^il + z)-^/^e-\ (6) 

where is the present age of the Universe. The Lorentz 
factor of the relevant string segment near the cusp is 
7 ~ 1/9. The duration of the cusp event as seen by a 
distant observer is 

r, ~ (1 + z)(at/2)7-3 ~ (aio/2)(l + z)-^/^e\ (7) 

One can expect that the observed duration of GRB is 
tgbb ^ Tc- This expectation will be justified by the 
hydrodynamical analysis below. 

The fluence, defined as the total energy per unit area 
of the detector, is 

S ^{l + z)(d£,^/dn)dl{z), (8) 

where = 3to{l+zY/'^[{l+zY/'^~l] is the luminosity 

distance. 

The rate of GRBs originating at cusps in the redshift 
interval dz and seen at an angle in the interval d6 is 
given by 

dNGRB ~ fs ■ ^edOil + z)-'iy{z)dV{z). (9) 

Here, ^{t) ^ ni{t)/Ti ^ 2a^^t^'* is the number of cusp 
events per unit spacetime volume, Ti ~ at/2 is the oscil- 
lation period of a loop, dV = 547rii^[(l + zY^^-lf{l + 
z)^^^/'^dz is the proper volume between redshifts z and 
z + dz, and we have used the relation dt^ = (1 -I- z)dt. 

Since different cusp events originate at different red- 
shifts and are seen at different angles, our model auto- 
matically gives a distribution of durations and fluences 
of GRBs. The angle 9 is related to the Lorentz factor of 
the relevant portion of the string as 9 ~ I/7, and from 
Eqs.(^,(||) we have 

7(z; S) ~ 7oa:85i//s:?^'[(A/rT^ - lfVT+~zf'^. 

(10) 

Here, 70 ~ 190, q;_8 = a/10~^, and the fluence S 
and the magnetic field Bq are expressed as 5 = S'_8 • 
10"® erg/cm^ and B = B_7 • 10"^ G. 

Very large values of 7 Jmax, which correspond (for 
a given redshift) to largest fluences, may not be seen 
at all because the radiation is emitted into a too nar- 
row solid angle and the observed rates of these events 



2 



are too small. The minimum value 7(2; Smin) is deter- 
mined by the smallest fluence that is observed, Smin ~ 
2 • 10~^ erg/cm^. Another limit on 7, which dominates 
at small z, follows from the condition of compactness ||] 
and is given by 7 > 100 (see below). 

The total rate of GRBs with fluence larger than S 
is obtained by integrating Eq.(^ over 6 from ^maxi^) 
to 7^^(z;5) and over z from to miTi[zm', zb], with Zm 
from jmax{zm) — l{zm] S). For relatively small fluences, 
S-s, < Sc = OM{j.rnaxiO)a-s/lo)^B'ir^, zb < Zrn and we 
obtain 

NGBBi> S) ^ ^ dViz){l + zfj-^iz; S) 

^ 3 -10^ SZl^^ yr-\ (11) 

Remarkably, this rate in our model does not depend on 
any string parameters and is determined (for a given 
value of S) almost entirely by the magnetic field Bq. The 
predicted slope NGB.Bi> S) oc S~^^^ is in a reasonable 
agreement with the observed one Nobs{> S) cx S~'^'^^ at 
relatively small fluences |l4[ . 

For large fluences S^s > Sc, integration of Eq.(||) gives 
Ngrb{> S) oc S~^/^. Observationally, the transition to 
this regime occurs at S^s ~ 10^ — 10"^. This can be 
accounted for if the cusp development is terminated by 
small-scale wiggles with fractional energy in the wiggles 
e ^ 10~'^aLgB't^j . Alternatively, if "fmax is determined 
by the back-reaction of the charge carriers, Eq.(||), then 
the regime (|l^) holds for larger S^s, and observed steep- 
ening of the distribution at large S can be due to the 
reduced efficiency of BATSE to detection of bursts with 
large 7. Indeed, large 7 results in a large Lorentz factor 
7CD of the emitting region (see below), and at •jcd ^ lO'^ 
photons start to escape from the BATSE range. 

The duration of a GRBs originating at redshift z and 
having fluence S is readily calculated as 

TGRB^200—f-^{l + z)-\Vrf^-l)-^s (12) 

O-S 

Estimated from Eq.(0), Tq^^ ^ 10^q;_8 s, while from 
Eq.(|l2|) using Smax ~ 1 • 10""* erg/cm^ and z ^ zb ^ 
4, one obtains Tq]:^ ^ 3algi?^Y ms. This range of 
TGRB agrees with observations for a_8 ^ 1 (i.e. 77 ^ 
2 • 10" GeV) and B_7 ~ 1. 

Let us now turn to the hydrodynamical phenomena in 
which the gamma radiation of the burst is actually gen- 
erated. The low-frequency electromagnetic pulse inter- 
acting with surrounding gas produces an ultrarelativistic 
beam of accelerated particles. This is the dominant chan- 
nel of energy loss by the pulse. The beam of high energy 
particles pushes the gas with the frozen magnetic field 
ahead of it, producing an external shock in surround- 
ing plasma and a reverse shock in the beam material, as 
in the case of "ordinary" fireball (for a review see [Q). 



The difference is that the beam propagates with a very 
large Lorentz factor 7f, ^ 7 where 7 is the Lorentz factor 
of the cusp. (The precise value of is not important 
for this discussion; it will be estimated in a subsequent 
publication |jl^.) Another difference is that the beam 
propagates in a very low-density gas. The beam can be 
regarded as a narrow shell of relativistic particles of width 
A ~ V^T^ in the observer's frame. 

The gamma radiation of the burst is produced as syn- 
chrotron radiation of electrons accelerated by external 
and reverse shocks. Naively, the duration of synchrotron 
radiation, i.e. tgrb, is determined by the thickness of 
the shell as tgrb ~ A. This is confirmed by a more de- 
tailed analysis, as follows. The reverse shock in our case 
is ultrarelativistic (l^Jl|. The neccessary condition for 
that, pf,//9 < 7^, is satisfied with a wide margin (here pb 
is the baryon density in the beam and p is the density of 
unperturbed gas). In this case, the shock dynamics and 
the GRB duration are determined by two hydrodynami- 
cal parameters [0. They are the thickness of the shell A 
and the Sedov length, defined as the distance travelled 
by the shell when the mass of the snow-ploughed gas be- 
comes comparable to the initial energy of the beam. The 
latter is given by Ised ^ {EisolpY^^- 

The reverse shock enters the shell and, as it propa- 
gates there, it strongly decelerates the shell. The syn- 
chrotron radiation occurs mainly in the shocked regions 
of the shell and of the external plasma. The surface sep- 
arating these two regions, the contact discontinuity (CD) 
surface, propagates with the same velocity as the shocked 
plasma, where the GRB radiation is produced. 

The synchrotron radiation ceases when the reverse 
shock reaches the inner boundary of the shell. This 
occurs at a distance i?A ^ ^sed^^^^ when the Lorentz 
factor of the CD surface is 7cd ~ (/sed/A)'^/* '-^ 
O-lB'^J^nZl'^il + ^)^/27^/^ where n_5 is the ambient 
baryon number density in units of 10~^ cm~^. Note that 
these values do not depend on the Lorentz factor of the 
beam 7f, and are determined by the cusp Lorentz factor 
7. The size of the synchrotron emitting region is of the 
order i?A, and the Lorentz factor of this region is equal to 
7ci5. The compactness condition requires Jcd ^ 10^, 
which yields 7 ^ 10^ used above. The duration of GRB 
is given by 

tgrb - Ra/2-i%o - (13) 

i.e. it is equal to the duration of the cusp event given by 
Eq. (0) . The energy that goes into synchrotron radiation 
is comparable to the energy of the electromagnetic pulse. 

We conclude with several remarks on our GRB model. 

(i) Magnetic fields in our scenario could be generated 
in young galaxies during the bright phase of their evolu- 
tion and then dispersed by galactic winds in the in- 
tergalactic space. One expects that at present the fields 
are concentrated in the filaments and sheets of the large- 
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scale structure lijjQ . With sheets of characteristic size 
L ~ (20 - 50)/i"^ Mpc and thickness D ~ S/i"^ Mpc, we 
have Jb ^ D/ L ^ 0.1. The average field strength can be 
estimated from the equipartition condition as Bq ~ 10~^ 
G ijl^ . We expect these values to be accurate within an 
order of magnitude. 

(ii) We emphasize that our model involves a number 
of simplifying assumptions. All loops at cosmic time t 
were assumed to have the same length I ^ at, while in 
reality there should be a distribution n{l,t). The evolu- 
tion law for B{z) and the assumption of Jb = const 
are also oversimplified. A more realistic model should 
also account for a spatial variation of B. As a result, the 
correlation between t and S suggested by Eq.(12) will 
tend to be washed out. Finally, we disregarded the pos- 
sibility of a nonzero cosmological constant which would 
introduce some quantitative changes in our estimates. 

(in) Apart from the reverse shock instability, small- 
scale wiggles on strings can naturally produce short-time 
variation in GRBs. These wiggles, acting like mini-cusps, 
produce a sequence of successive fireballs, which are the 
usual explanation of GRB short-time structure [Q. 

(iv) Cusps reappear on a loop, producing nearly identi- 
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yr. 



In 



cal GRBs with a period T; - at/2 - 50(H-z) 
a more realistic model, some fraction of loops would have 
lengths smaller than at and shorter recurrence periods. 

(v) Our model predicts that GRBs should be accom- 
panied by strong bursts of gravitational radiation. The 
angular distribution of the gravitational wave energy 
around the direction of the cusp is |^ dSg/dfl ^ 
and the dimensionless amplitude of a burst of duration r 
originating at redshift z can be estimated as 



kg{G^lf{T/ato)'^'{l + z)-'/'[il + _ 1] 



(14) 



or h 



lO-^^a^J^z- 



^{r/lsy/^ for z < 1. Here, we have 
used the relation Fg - h^/Gr^ - (1 + z)(d£g/dn)dl^{z) 
for the gravitational wave flux and Eq.(|j) for the burst 
duration t. These gravitational wave bursts are much 
stronger than expected from more conventional sources 
and should be detectable by the planned LIGO, VIRGO 
and LISA detectors. 

(vi) GRBs have been suggested as possible sources 
of the observed ultrahigh-energy cosmic rays (UHECR) 
[^,|l|. This idea encounters two difhculties. If GRBs 
are distributed uniformly in the universe, UHECR have 
a classical Greisen-Zatsepin-Kuzmin (GZK) cutoff, ab- 
sent in the observations. The acceleration by an ultra- 
relativistic shock is possible only in the one-loop regime 
(i.e. due to a single reflection from the shock) For a 
standard GRB with a Lorentz factor ^gh ^ 300 it results 
in the maximum energy E„iax ^ llh'^p ^ ^^^^ ^^'^ 
too low for UHECR. 

Our model can resolve both of these difhculties, assum- 
ing that gammamax is determined by the current back- 



reaction, Eq.(||). 

If the magnetic field in the Local Supercluster (LS) is 
considerably stronger than outside, then the cusps in LS 
are more powerful and the GZK cutoff is less pronounced. 

Near-cusp segments with large Lorentz factors pro- 
duce hydrodynamical flows with large Lorentz factors, 
e.g. -"y ~ 2 • 10"' corresponds to ^cd ^ 3 • 10^ and 
Emax ~ IcD^p ~ 1 ' 10^" ■ Protons with such en- 
ergies are deflected in the magnetic field of LS and can 
be observed, while protons with much higher energies 
caused by near-cusp segments with 7 ^ 10^ are propa- 
gating rectilinearly and generally are not seen. Further 
details of this scenario will be given elsewhere ||l5[] . 
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